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CONVERGENCE OF THE FULL COMPRESSIBLE 
NAVIER-STOKES-MAXWELL SYSTEM TO THE 
INCOMPRESSIBLE MAGNETOHYDRODYNAMIC EQUATIONS 
IN A BOUNDED DOMAIN 

JISHAN FAN, FUCAI LI* *, AND GEN NAKAMURA 

Abstract. In this paper we establish the uniform estimates of strong solu¬ 
tions with respect to the Mach number and the dielectric constant to the full 
compressible Navier-Stokes-Maxwell system in a bounded domain. Based on 
these uniform estimates, we obtain the convergence of the full compressible 
Navier-Stokes-Maxwell system to the incompressible magnetohydrodynamic 
equations for well-prepared data. 


1. Introduction 

Ill this paper we consider the singular limit of the following full compressible 


Navier-Stokes-Maxwell system in a bounded domain C (0): 

+ div (pu) = 0, (1.1) 

dt{pu) + div {pu 0 m) H— 2 ^P — pAw — (A + p)VdivM 

= {E + uxh)xb, (1.2) 

dt{pe) + div (pue) + pdiv u — div (kVT) 

= ei{2p,\D(u)\‘^ + A(divM)^ + {E + u x b)^), (1.3) 

e2dtE — Totb + E + UX 6 = 0, (1.4) 

9(6 + rot £’ = 0, div6 = 0, (1.5) 


where the unknowns pyU^PyCyE, E, and 6 stand for the density, velocity, pressure, 
internal energy, temperature, electric field, and magnetic field, respectively. The 
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physical constants ^ and A are the shear viscosity and bulk viscosity of the flow and 
satisfy /r > 0 and A + |/r > 0. k > 0 is the heat conductivity. ei > 0 is the (scaled) 
Mach number, and e 2 > 0 is the (scaled) dielectric constant. D{u) := i(Vu+Vu"'"), 
where denotes the transpose of the matrix Vit. 

In [8l|9], Kawashima and Shizuta established the global existence of smooth 
solutions for small data m and studied its zero dielectric constant limit e 2 —>■ 0 in 
the whole space R^. Recently, Jiang and Li [Q studied the zero dielectric constant 
limit €2 —>■ 0 to the system (fLT])-(IT3D and obtained the convergence of the system 
(fm)-(fT31) to the full compressible magnetohydrodynamic equations in , see also 
[7] on the similar results to the invisid case of (ll.II) - H1.5D . In [TU], Li and Mu study 
the low Mach number limit ei —>■ 0 to the system (frT])-(IT31) and obtained the 
convergence of the system (HHl-dESl) to the incompressible Navier-Stokes-Maxwell 
system in the torus T^. 

It should be pointed out that no boundary effect is considered in the references 
mentioned above. The purpose of this paper is to invistigate the singular limit 
Cl, £2 —>■ 0 to the system (fT3])-(IT31) in a bounded domain. For simplicity, we shall 
take £1 = £2 = e and consider the case that the fluid is a polytropic ideal gas, that 
is 

e := CuT, p := TZpT (1.6) 

with Cv > 0 and TZ being the specific heat at constant volume and the generic gas 
constant, respectively. 

To state the main result of this paper, we denote the density and temperature 


variations by cr'^ and 

:= l + ea% := l + £dL (1.7) 

Then we can rewrite the system (frT])-(IT31) as follows: 

+ div (fj'^it'^) + -div = 0, (1.8) 

£ 

p^{dtu^ + • Vu^) + y(Va'^ + - (A + /i)Vdiv 

= {E^ + X b'^) X b\ (1.9) 

Cvp^{dtd^ + ■ V9'^) + Tl{p’^9'^ + (T'^)divit'^ + —divu*^ 

€ 

= kA 9^ + e[2p\D{u^)\^ + A(divu'^)^ + + u^x E)^], (1.10) 

edtE'^ — rot E + E'^ + u'^ x E = 0, (l-H) 
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dtb^ + rot = 0, div E = 0. (1.12) 

Here we have added the superscript e on the unknowns (cr, u, 0, E, b) to emphasise 
the dependence of e. The system (I1.8I) - (I1.12I) are supplemented with the following 
initial and boundary conditions: 


in fl, (1.13) 

Qoe 

u'^ ■ n = 0, rotu*^ x n = 0, = 0, E'^ X n = 0, E ■ n = 0 on 911, (1-14) 

an 

where n is the unit outer normal vector to the smooth boundary dfl. 

Formally, if we let e —>■ 0 in dm and (HID, then we obtain that div —>• 
0, —>■ 0, and Vcr'^ —>■ 0. Letting e = 0 in (11.111) gives E'^ = lotE — u'^ x E. 

Pulling it into (11.121) and taking the limit e —>■ 0 we obtain the following limit 
system (suppose that the limits —>■ {v,B) exist): 

vt + v- Vv + Vtt — = rot B x B, 

Bt+Tot{Bxv)-AB = 0, (1.15) 

divu = 0, div 5 = 0. 

We shall give a rigorous proof the the above formal analysis below. 

Denote 

M'^(t) := sup |||(cr^u^r,ye^;^5'^)(•,s)||^^2 + ||9t(cr% r, Ve-E", 5'^)(-, s)||^i 

0<s<i ^ 


+ e||9;(cr",'u",6»'^)(-,s)||i2 

rt 


.} 


1 + ecr'^(-, s) 

+ I (||(uS0^)||?,3 + + \\edUa^u^,e^)rm 


(1.16) 


First, we have 


Theorem 1.1. Let 12 C 6e a simply connected, bounded domain with smooth 
boundary 912 and 0 < e < 1. Suppose that the initial data ((t% 0% 5% 6g) satisfy 

the following regularity conditions: 


o<0^, 


0 < —— < 1 + eo-Q < Kq, 
Eo 


To,Uo,0^o,E^,Eq)\\h2 + ||9t(CT%M%6»%5%6%(-,0)||ifi 


(1.17) 


+ e||9,^(cr%M%5)(-,0)||i2 < Ki 


(1.18) 
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for some positive constants Kq > 1 and Ki independent o/e > 0. Then there exist 
a small time T > 0 independent of e > 0 and a unique strong solution (cr, u, 0, E, b) 
to the initial boundary value problem (I1.8D - (I1.14I) such that 

M^{f) < K (1.19) 

for some positive eonstant K independent of e> 0. 

Remark 1 . 1 . In the assumption (11.181) . crj(-,0) is indeed defined by —div (ctqUq)+ 
■idivug through the density equation and the other quantities are defined by an 
analogous way. 

Based on the uniform estimates of the solutions, we can prove the following 
convergence result by applying the Arzela-Ascolis theorem in a standard way. 

Theorem 1 . 2 . Let , u'^, 6'^, E'^, E) be the solution of the problem (I1.8I) - (I1.14I) with 
initial data ((Tq, rtg, 0g, Ag, &g) satisfying the conditions in Theorem ll.il Assume 
further that the initial data {aQ,UQ,0Q, Eg^bl) satisfy that 

(ecTg, Ug, e0g, &g) —>■ (O, Vq , 0, Bq) sEongly in for any 0<s<2 as e—>-0, 

Eq —>• rot Bq — Vq X Bq strongly in iJ® for any 0<s<l as e—>-0. 

Then (etr®, m®, E) —>■ (0, v, 0, B) strongly in L°°{0, T ; H^) and E'^ ^ B — v x B 

strongly in L°°{0,T; L^) as e 0, where (y,B) satisfies (11.151) with the following 
initial and boundary conditions: 

{ V ■ n = B ■ n = 0, rot v x n = rot B x n = 0 on dft x (0, Tl, 

^ ^ (1.20) 

{v,B){-,0) = {vo,Bo){-) in C 

The remainder of this paper is devoted to the proof of Theorem 11.11 which will 
be given in next section. 


2. Proof of Theorem O 

In this section we shall prove Theorem 11.11 by combining the ideas developed 
in [HEiiiiin]. First, by taking the similar arguments to that mm, we know that 
in order to prove (11.191) . it suffices to show the following inequality 

M®(<) < Co(M®(0)) exp[t3C'(A/®(t))] (2.1) 

for Vt G [0,T] and some given positive nondecreasing continuous functions C'o(-) 
and C'(-). 
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Below we shall omit the spatial domain in the integrals and drop the super¬ 
script “e” of 0*^, etc. for the sake of simplicity; moreover, we write M^{t) 

and M^(0) as M and Mq, respectively. Since the physical constants k, Cy, and TZ do 
not bring any essential difficulties in our arguments, we shall take k = Cy = 7^ = 1. 
We will also use the following two inequalities: 

||M||ff'»(n) < C'dldivullffo-qn) + ||rotM||/i-<.-i(o) -I- HullLi-o-qn) + \\u ■ n||^«-i/2(9Q)), 


( 2 . 2 ) 


l|M||ff'>(n) < C'(l|divM||^/«-i(n) -f ||rotM||^/«-i(n) -f ||m||^«-i(o) -h \\u x n||//«-i/ 2 (gn)). 


(2.3) 


for any u G with s > 1, which were obtained in [2] and [12] respectively. 

Because the local existence for the problem (I1.8I) - (I1.14I) with fixed e > 0 is 
essential similar to that in m, we only need to prove m- We will use the 
methods developed in mM- 

First, by the same calculations as that in [3], we get 


-(•:0 +Ilp(-:0llff2 <Co{Mo)exp{CVtM) 

p L~ 


(2.4) 


(2.5) 


Now we use the same method as that in |4] to prove some a priori estimates on 


(if, 6). 


Testing (11.111) and (11.121) by E and b, respectively, and summing up the results, 
we see that 




Integrating the above inequality over (0,t), we find that 



( 2 . 6 ) 


Using (11.141) and the formula 


{u X b) X n = {b ■ n)u — {n ■ u)b = 0 on 917, 


(2.7) 


we infer that 


rot b X n = 0 on 9f7. 


(2.8) 
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Taking rot to (jl.llj) and (jl.l2|l . testing the results by rotE and rot6, respec¬ 
tively, summing up the results, and using ()2.8I1 and integration by parts, we have 


J + J |roti?pfia; 

= — J rot {u X b) ■ rot Edx 

<\j \T:otE\^dx + C\\u\W-2\\ioih\\l.2 

< - /" |rot iilpfia; -I- C(M). 


Integrating the above inequality over (0,t), we have 

y (e|rot Sp -I- |rot -I- J J \iot E\^dxds < Co{Mq) + tC{M). (2.9) 
Taking div to (11.111) and testing the result by div E, we infer that 
J + J (divi?)^(ia; = J div (6 x u)dWEdx 

<^J {div Efdx + C{M). 

Integrating the above inequality over (0,t), we deduce that 

e J{divE)^dx + J J {div E)^dxds < Cq{Mq) + tC{M). (2-10) 

Taking dt to (11.111) and (11.121) . testing the results by Et and bt, respectively, 
summing up the results, we get 
1 d 


2 dt 


\e\Et\'^ + \bt\^)dx + J \Et\^dx = J dt{b X u)dtEdx 


< ly \Et\^dx + C{M). 

Integrating the above inequality over (0,t), we get 

J{e\Et\^ + \bt\^)dx + y y \Et\^dxds < Co(Mo) + tC{M). (2.11) 


(11.121) and (12.81) give the boundary condition 

rot '^E X n = 0 on dfl. 


( 2 . 12 ) 


Taking rot^ to (11.111) and (I1.12L testing the results by rot^i? and rot ^6, respec¬ 
tively, summing up the results, we derive that 
1 d 


2 dt 


J{e\rot'^E\'^+ \rotH\^)dx + J \rot'^E\^dx 
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= J Tot^{b X u)rot‘^Edx 

-\j + 

Integrating the above inequality over (0,t), we have 

y(e|rot^ii^P + |rot^&P)dx + J J liot"^E\‘^dxds < Co{Mo) + tC{M). (2.13) 

Taking Vdiv to (ll.llll . testing the result by Vdivi?, we have 

J J l^divii^pda; = J Vdiv {b x u)-Vdiv Edx 

-\j l^divi;|2dx + C'||6||^2||u||^2 

<^J \VdivEfdx + C{M). 

Integrating the above inequality over (0,t), we obtain 

e J |Vdivi?pda: + J J \V div E\‘^dxds < Co{Mo) + tC{M). (2-14) 

Taking 9trot to (11.1111 and (I1.12L testing the results by dtrotE and dtTOtb, 
respectively, summing up the results, and using (12.811 . we obtain 

y(ekoti^tP + |rot6tp)da: + y|rot£'tpda; 

= J Tot {bt X u + b X ut)rot Etdx 

< C{\\bt\\m\\u\\H^ + \\b\\H2\\ut\\m)\\iotEt\\L2 

— 2J + C{M). 

Integrating the above inequality over (0,t), we obtain 

y(e|rot + |rot6t p)(ia: + J J {rot Et\^dxds < Co{Mq) + tC{M). (2.15) 

Applying 9tdiv to (jl.lljl . testing the result by divE^j, we have 
e d 


2 dt 


y (divEt)^da; + J{div Et)^dx 
= y div {bt X u + b X ut)div Etdx 


< C'(||6t||/i-i||u||jy2 + Iluilli/i||6||ij2)||divEt||i2 

< ^ J {div Etfdx + C{M). 
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Integrating the above inequality over (0,t), we have 


e J {div Et)‘^dx + J j {ddw Etfdxds < Cq{Mq) tC{M). (2.16) 


Now we use the method in [3] to prove some a priori estimates on (a,u,d). 
Testing (lEHD, (HH) and (11.101) by a, u and 9, respectively, summing up the results, 
we obtain 


l_d 
2 dt 


J {<7^ + pu^ + p9^)dx + J {p\Vu\'^ + (A + ^)(div + \V9\‘^)dx 
= J div u — pO^'^ dx + e J 6{2p\D(u)\'^ + X{div u)"^ + {E + u x b)'^)dx 


<\\ynU^{Mh + \\p\\L^mh) 

< + C||Vu|U~(||Vu|U. + \\E\\l, + ||u||io„|16||i4)||d|U. 
<\\Vu\\l^C{M) < \\u\\h^C{M). 


Integrating the above inequality over (0,t), we obtain 


(cT^ + + p9‘^)dx + J (|Vu|^ + \Ve\^)dxds < Co{Mq) + ViC{M). (2.17) 


Applying dt to (11.81) . (11.91) and (11.101) . we see that 

dtt + -divut = -div (cru)t, 
e 

p{uu + u ■ Vut) + -(VcTt + VOt) - p/Xut - (A + p)Vdiv Ut 
e 

= -ptut - {pu)tVu - V(cr9)t + [(E + u X b) X b]t, 


(2.18) 


(2.19) 


p{9tt + u ■ V9t) + -div Ut - A9t = -ptdt - {pu)t ■ S/O - {{p6 + a)divu)t 


+ e{2p\D{u)\‘^ + A(divu)^ + {E + u x b)'^)t. 


( 2 . 20 ) 


Testing (I2.18p . (12.191) and (12.201) by <7t,ut and d*, respectively, summing up the 
results, we reach 

2dt J d-put + p9^)dx + J {p\S/ut\'^ + {X + p){div ut)'^ + \S/0t\'^)dx 
= J {au)tyatdx - J [ptUt + [pu)tVu + V{a9)t]utdx 
+ J {{E + uxb) X b)tUtdx - J [pt9t + {pu)tV9 + {{p9 + cr)div u)t]9tdx 
d~ ^ J (2p|71(u)|^ + A(divtt)^ + {E + u X b)‘^)t9tdx 

< C{M) + J {{E + uxb) X b)tUtdx d- ^ J {{E + ux b)‘^)t9tdx 
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<C{M) + \\Et\\L^C{M). 

Integrating the above inequality over (0,t), we have 
j {a^ + pu^ + p9^)dx + J J {\Vut\^+ \V 6 t\^)dxds < Co{Mo) + ViC{M). (2.21) 
Testing (12.191) by —Vdivu in x (0,t)), we find that 


2 

/r + A 


^ II Vdivu(-, t )|||,2 - J J {Vat + 'V9t) ■ Vdiv udxds 

||Vdivuo||i2 + J J p{utt + u ■ Vut)Vdiv udxds 
+ J J {ptUt + {pu)tyu + {a9)t) ■ Vdivudxds 
— J J {{E + u X b) X 6)tVdiv udxds 


p X. 


Vdiv MqIIls + di + /2 + /a- 


( 2 . 22 ) 


We bound /i ,/2 and I 3 as follows. 


h<C{M) [ \\utt\\L^ds + tC{M)<ViC{M), 
Jo 


h <tC{M), 


h <C{M) [ \\Et\\H^ds + tC{M) < ViC{M). 

Jo 

Applying V to (11.81) and (II.IOL testing the results by Vert and V0t in L^(r2 x 
(0,t)), respectively, we derive 


^dxds 


VtTtVdiv udxds 


= — J J Vd\v {au)Vatdxds <tC{M), 


(2.23) 


and 


- J \A9f dx + J J p|V6)t|^da;ds H —J J V9t'X/div udxds 

= iy \A9o\'^dx + J J eV[2p\D{u)\^ + X{d\vu)^ + {E + u X b)^]W9tdxds 
— J J V[(p0 + (T)div M]V0tda;ds — y J [Vp9t + V{pu ■ V9)]V9tdxds 


< Co{Mo) + ViCiM). 


(2.24) 
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Summing up (I2.22p . (I2.23|) and (I2.24|l . we arrive at 

J (iVdivMp + \A9\^)dx + J J (|Vatp + \Vet\‘^)dxds < Co(Mo) exp(v^C'(M)). 

° (2.25) 

Testing (I2.18|) . (12.191) and (j2.20|) by —Aat, — Vdivit* and —A6t, respectively, we 
derive 


\Vcrt\'^dx + ^ I I V dW UtV (Jtdxds 


— J |Vo-t(0)p(ia; + J J div {cTtU + <JUt) ■ Aatdxds 
:^J\Vatm^dx + h. (2.26) 


We bound as follows. 

rt 


li = J J uVatAatdx — J J V{atdiv u + UtVcr + crdiv Ut)ycrtdxds 
= — J J diuVatdiCTtdxds + J J — div u\Vat\^dxds 

— J J {atdiv u + Ut\7a + adiv Ut)\7atdxds 

<C{M) f WVuU’^ds + C{M) [ \\u\\H^ds + C{M) [ \\ut\\H^ds 
Jo Jo Jo 


<VtC{M). 


And 


i J p{divut)^dx + {X + 2^) J J |Vdivutpda;c?s 
~~ J J Vdiv Mt(V(Tt + V0t)(ixds 

= ^ J po{divut{0))^dx + J J[ptUt + {pu)tX/u + X/{a9)t]X/div utdxds 
+ J J div (i? X b)tdiv utdxds — J J [{u x b) x b]tydiv utdxds 
+ J J{eV a ■ utt + u ■ y Ut)div Utdxds — J J VuidjUtdiv Utdxds 

<Co(A/o) + VtC(A/). (2.27) 

+ J J\A9t\^dxds + - J J y div utV 9tdxds 


And 

5 
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= ^ y Po|V0t(O)p(ia; - y J eWa9ttV9tdxds - J J 'Y'y {pUi)di9tS/9tdxds 

+ J J A9[pt9t + {pu)t'V9 + {{p9 + a) ■ div u)t]dxds 

~^J y A9t{2p\D{u)\‘^ + X{div u)'^ + {E + u X b)'^)tdxds 
<Co{Mo) + ViC{M). (2.28) 

Summing up ()2.26p . (I2.27P and ()2.28p . we arrive at 

y i\Vat\y {divutf + \V9tf)dx + j{\VdWutf + iA9tf)dxds 


< C'o(Mo) exp(v^C'(M)). 


(2.29) 


Now, testing di\7 (II.8p by di\7a + diV9 and the same calculations as those in [5] 
to obtain 

i y \diX/a\'^dx + J diVa ■ diV9dx + - J J diX/divu{diX/a + diX/9)dxds 

<Co{Mo) + ViC{M). (2.30) 

Testing di (11.91) by ^^Vdivu in T^(r2 x (0,t)) and the same calculations as those 
in [3] to obtain 

iy y l^iVdivMpdx — - y y di'Vdivu{di'Va + di'V9)dxds<tC{M). (2.31) 
(|2.30l) . (12.311) and (12.251) give 

y iV^crpda; + J J div u\‘^dxds < Co{Mo) exp{\/tC{M)). (2.32) 

Applying rot to (11.91) and denoting the vorticity ui := lotu, we see that 

pujt + pu-Vuj - pAuj ={djpuit - dipujt) + {dj{puk)dkUi - di{puk)dkUj) 

+ rot [{E + u X b) X b]. (2.33) 

We test (12.331) by Aw in i^(n x (0,t)) to get 

y p\iot uj\'^dx + p J J \Aijj\^dxds < Co{Mtj) + \/tC{M). (2.34) 

Similarly, we apply dt to ()2.33l) and test the resulting equations by ujt in L^(r2 x 
(0,t)) to deduce that 

y p\uJt\^dx + p J J \rotLUt\‘^dxds < Co{Mq) + \/tC{M). 


(2.35) 
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By the same calculations as that in [3], we have 

l|A 6 '||L 2 (o,t./fi) < C'o(Mo)exp(t3C'(M)). (2.36) 

It follows from (fTTll . (fTT^ . (1^ . (IZ^ . (I2l0ll . (imi) . (I2T^ . (l214ll . (l2J5ll 
and (12.161) that 

e||'E'tt||L 2 (o,t;L 2 ) < ||rot6t — Et — {u x b)t\\L^(^o^t;L^) < Co{Mq) + tC{M). (2.37) 

V~4bttit)\\L^ < V~4^otEt{-ML^ < C'o(Mo) + tC(M). (2.38) 

Finally, we need to estimate tott, iutt and eOtt to close the energy estimates. 
Testing g^^ dl.SI) . cl? (ll.9|) and (II. 101) by e^au, e^Utt and respectively, then 

by the same calculations as that in [3], we conclude that 

4{<^tt,utt,0tt){t)\\L^ +e\\iutt,0tt)\\L^io,t;H^) < Cq(Mq) exp(t^C'(M)) 

and thus (11.191) hold true. 

This completes the proof the proof of Theorem ll.il 

□ 
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